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Very high energy (VHE, energy E > 100 GeV) y-rays from cosmological sources are attenuated due to the interaction with photons 
of the extragalactic background light (EBL) in the ultraviolet to infrared wavelength band. The EBL, thus, leaves an imprint on the 
observed energy spectra of these objects. In the last four years, the number of extragalactic VHE sources discovered with imaging 
atmospheric Cherenkov telescopes (lACTs), such as MAGIC, H.E.S.S., and VERITAS, has doubled. Furthermore, the measurements 
of the Fermi satellite brought new insights into the intrinsic spectra of the sources at GeV energies. In this paper, upper limits on the 
EBL intensity are derived by considering the most extensive VHE source sample ever used in this context. This is accomplished by 
constructing a large number of generic EBL shapes and combining spectral informations from Fermi and lACTs together with minimal 
assumptions about the source physics at high and very high y-ray energies. The evolution of the EBL with redshift is accounted for 
and the possibility of the formation of an electromagnetic cascade and the implications on the upper limits are explored. The EBL 
density at z = is constrained over a broad wavelength range between 0.4 and 100 [xm. At optical wavelengths, the EBL density is 
constrained below 24 nW m"- sr" ' and below 5 nW m"^ sr" ' between 8 [rm and 3 1 ^m. 
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1. Introduction 

The extragalactic background light (EBL) is the diffuse and 
isotropic radiati on field from ultraviolet to far infrared wave- 
lengths (see e.g. l Hauser&Dwe3l200ll: lKashlinskvll2005 . for a 
review). It originates from the starlight integrated over all epochs 
and starlight emission reprocessed by interstellar dust. These 
two distinct components lead to two maxima in the spectral en- 
ergy distribution (SED) of the EBL, the first at ~ 1 [xm (starlight) 
and the second at ~ 100 |,im (dust). Narrow spectral features like 
absorption lines are smeared out in the integration over redshift 
leading to a smooth shape of the EBL at z = 0. Further con- 
tributions may come from diffuse emission from galaxy clus- 
ters (Chelouche et al. 2007), unresolved act ive galactic nuclei 
("Matute et al. 2006), Population III stars (e.g. lSantos et al.l2002t 
Raue et al. 2009), or exotic sources like dark m atter powered 
stars in the early universe (iMaurer et al 1 12011' ). The SED of 
the EBL at z = comprises information about the star and 
galaxy formation rates and the dust content of galaxies. Direct 
measurements of the EBL are, however, impeded, especially in 
the infrared, by for eground emission such as the zodiacal light 
dHauser et al. I [l998l) . Therefore, upper and lower limits are of- 
ten the only available information about the EBL density. Lower 
limits are derived from integrated galaxy num ber counts e.g. by 
the Hubble Space Telescope in the optical (Madau & Pozzettil 
2000) and the Spitzer telescope in the infrared (Fazio et aL 
20041) . Several authors have modeled the EBL in the past 



(e.g.lPrim ack et a l.llIOOllStecker e t al."2006l: 'Franc eschini et al.1 



l2008H Kneiske & bolell2010l: lD'ommguez et al. 201 ll). Although 



the approaches forecast different EBL densities at z = 0, the 
most recent models more or less agree on the overall EBL shape. 
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The observations of very high energy (energy E > 100 GeV; 
VHE) y-rays from extragalactic sources with imaging atmo- 
spheric Cherenkov telescopes (lACTs) has opened a new win- 
dow to constrain the EBL density. Most of the extragalactic y- 
ray sources are active galactic nuclei (AGN), especially blazars 
(see e.g. Urrv & Padovani 1995). The y-rays from these ob- 
jects are attenuate d by the pair production mechanism: yvHE + 
Tebl ^e++e- (N ikishovll962H.Tellevlll966HGould & Schreded 
1967). If assumptions are made about the properties of the in- 
trinsic blazar spectrum, a comparison with the observed spec- 
trum allows to pla ce upper limits on the EBL intensity (e.g. 
IStecker et al.lll992l) . In this context, th e spectra of Markarian 
(Mkn) 501 during an extraordinary flare (^Aharonian et al.''1999) 
and of the distant blazar H 1426-1-482 ( Aharonian et al. 20031) 
resulted in the first constraints of the EBL density from mid to 
far infrared (MIR and FIR) wavelengths. With the new genera- 
tion of lACTs, limits were derived from the s pectra of lES 1 101- 
232 and H 2356-309 dAharonian et sH IOOd) and lES 0229-H200 
(Aharo nian et al.ll2007a ) in the near infrared (NIR), and in the 
optica l from t he MAGIC observation of 3C279 (Albe rt et all 
20081) . lMazin"& Raue (2007, henceforth MR07) use a sample of 
all at that time known blazars and test a large number of differ- 
ent EBL shapes to derive robust constraints over a large wave- 
length range. The authors exclude EBL densities that produce 
VHE spectra, characterized by dN/dE oc E^^, with F < Fiimit 
(being Fiimit =1.5 for realistic and Fumit - 2/3 for extreme sce- 
narios) or an exponential pile up at highest energies. 

With the advent of the large area tel escope (LAT) on board 
the Fermi satellite (lAtwood et al.ll20()9l) and its unprecedented 
sensitivity at high energies (100 MeV < £ < 100 GeV; 
HE), further possibilities arose to confine the EBL density. 
Bounds can be derived either by considering solely Fermi- 
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LAT observations of AGN and gamma-ray bursts (lAbdo et alJ 
2010b; Raue 2010") or by combining HE with VHE spectra (e.g. 
.Georganopoulos et al. 2010; Orretal. 2011). It has also been 
proposed that the Ferm/-LAT can, in principle, measure the EBL 
photo ns upscattered by electrons in lobes of radio galaxies di- 
rectly (iGeorganopoulos et al.ll2008l) . Attenuation limits can also 
be estimated by modeling the entire blazar SED in order to 
forecast the intrinsic VH E emission ( Krawczynski et al. 200^ 
iMankuzhivil et al.ll2010l) . 

In this paper, results from the recently published Fermi 
two year catalog (iThe Fermi-LAT collaboration! 1201 li hence- 
forth 2FGL) together with a comprehensive VHE spectra sample 
are used to place upper limits on the EBL density. This approach 
relies on minimal assumptions about the intrinsic spectra. The 
VHE sample is composed of spectra measured with different in- 
struments, thereby ensuring that the results are not influenced 
by the possible systematic bias of an individual instrument or 
observation. 

The article is organized as follows. In Section |2] the calcu- 
lation of the attenuation is presented in order to correct the ob- 
served spectra for absorption. The resulting intrinsic spectra are 
subsequently described with analytical functions. Section|3]out- 
lines in detail the different approaches to constrain the EBL be- 
fore the selection of VHE spectra is addressed in Section |4] The 
combination of VHE and HE spectra of variable sources will 
also be discussed. The results are presented in Section |5] be- 
fore concluding in Section |6] Throughout this work a standard 
ACDM cosmology is assumed with 0„, - 0.3, = 0.7, and 
h = 0.72. 



2. Intrinsic VHE gamma-ray spectra 

The intrinsic energy spectrum dA^int/d/? of a source at redshift 
Zo at the measured energy E differs from the observed spectrum 
dA^obs/d£ due to the interaction of source photons with the pho- 
tons of the EBL which is most commonly expressed as 



d£ 



d£ 



:exp[-Ty(£',zo)]- 



(1) 



The strength of the attenuation at energy E is given by the optical 
depth Ty-. a threefold integral over the distance {, the cosine fi of 
the angle between the photon momenta, and the energy e of the 
EBL photons (e.g. .Dwek & Ki'ennrich.2005i) . 



Ty{E, Zo) 



ZO +1 



J d€{z) j dju^ j d£'nEBL{€\z)(ryy(E',£',fi). (2) 



-1 



The primed values correspond to the redshifted energies and 
nEBL(f',z) denotes the comoving EBL photon number density. 
The threshold energy for pair production is given by e^^^^ - 
ethi(E',fi) with E' = £(1 + z). The cross section for pair pro- 
duction, (Tyy, is strongly peaked at a wavelength (e.g. , Guv et al.> 
l2000h 



he / E \ 



(3) 



and, therefore, VHE y-rays predominantly interact with EBL 
photons from optical to FIR wavelengths. 

The comoving EBL photon density is described here by 
splines constructed from a grid in (/i,v/v)-plane (see MR07 for 




Fig. 1. Upper panel: Grid in wavelength versus the energy den- 
sity of the EBL used to construct the EBL shapes for testing 
(red bullets). Also shown are the minimum and maximum shape 
tested (solid lines) and the same for the grid of MR07 (blue trian- 
gles; dashed lines). Lower panel: Minimum and maximum EBL 
shape tested versus EBL limits a nd measurements (light gray 
symbols, see lRaue & Mazinll201 iL and references therein). 



further details). This ensures independence of EBL model as- 
sumptions and allows for a great variety of EBL shapes to be 
tested. Furthermore, the usage of splines drastically reduces the 
effort to compute the complete threefold integral of Eq.|2]numer- 
ically as shown in MR07. 

Each spline is defined by the choice of knot points and 
weights from the grid in the (/l,v/y)-plane. The grid is bound 
by a minimum and a maximum shape, shown in Figure [1] The 
setup of grid points is taken from MR07. The minimum shape 
tested is set to reproduce the lower li mits from the galaxy num- 
ber counts from Spitzer (iFazio et al.l [2004) while the maximum 
shape roughly follows the upper limits derived from measure- 
ments. To reduce the computational costs, the extreme cases con- 
sidered by MR07 of the EBL density in the optical and near in- 
frared (NIR) are not tested. Moreover, with current VHE spectra 
the EBL intensity is only testable up to a wavelength /I » 100 [xm 
so no additional grid points beyond this wavelength are used. In 
total this range of knots and weights allows for 1,920,000 dif- 
ferent EBL shapes. A much smaller spacing of the grid points 
is not meaningful as small structures are smeared out in the cal- 
culation of Ty and the EBL can be understood as a superposi- 
tion of black bodies that are not arbitrarily narrow in wavelength 
(MR07, Raue 2007). 

In most previous studies, no EBL evolution with redshift 
is assumed when computing EBL upper limits using VHE y- 
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ray observations. Neglecting the evolution leads to an overesti- 
mation of the optical depth between 10% (z = 0.2) and 35 % 
(z - 0.5) (Raue & Mazin 2008) and, consequently, too rigid up- 
per limits (see Appendix lAli. In this study, the z evolution is ac- 
count ed for by a phenomenological ansatz (e.g. iRaue & MazinI 
l2008h : the effective cosmological photon number density scal- 
ing is changed from «ebl (1 -h z)-' to hebl (1 + z)^^-^"°- For 
a value of /evo = L2 a good agreement is found between this 
simplified approac h and complete EBL model calculations for 
redshifts z < 0.7 dRaue & MazinI l2008h . Including the redshift 
evolution of the EBL in general decreases the attenuation com- 
pared to the no-evolution case and, therefore, weaker EBL limits 
are expected 

The intrinsic y-ray spectrum for a given EBL shape and mea- 
sured y-ray spectrum is reconstructed by solving Eq.[T]for dA^im/ 
dE. For a spectrum with « energy bins the relation reads 




where the energy of the logarithmic bin center is denoted by 
Ei. A systematic error is introduced by using Ty calculated for 
the energy at the bin center since, on the one hand, the attenua- 
tion can change dramatically within relatively wide energy bins 
and, on the other hand, the mean attenuation actuall y depends on 
the i ntrinsic spectral shape in the energy bin ( Stec ker & Scullvl 
12009 ). The introduced error is studied by comparing Ty with 
an averaged value of the optical depth over the highest energy 
bin for the spectra that are attenuated most. These spectra are 
described with an analytical function f{E) (a power or broken 
power law, cf. Table |2]i and the averaged optical depth (Ty) is 
found to be 

/ Ty(E,z)f(E)dE 

(Ty) = ^-^r ■ (5) 

/ f(E)dE 

AE 

The results are summarized in Table [T] The ratios {Ty}/Ty are 
close to, but always smaller than, one and the optical depth 
is overestimated by < 5 %. Thus, the simplified approach adds 
marginally to the uncertainties of the upper limits. 

Table 1. Comparison between the optical depth at the logarith- 
mic bin center of the highest energy bin and the averaged value 
over the bin width. 



Source 


Minimum EBL siiape tested 

Ty (Ty) {ry)/Ty 


Maximum EBL shape tested 

Ty (Ty) (Ty)/Ty 


3C279 


3.48 


3.34 


0.96 


18.33 


17.61 


0.96 


H 1426+428 


2.54 


2.53 


0.99 


12.61 


12.48 


0.99 


lES 1101-232 


2.69 


2.68 


1.00 


13.62 


13.57 


1.00 


MknSOl 


3.27 


3.21 


0.98 


11.86 


11.67 


0.98 



Notes. The HEGRA spectrum is used here for Mkn501, see Table[3] 
for the references. 

The intrinsic spectra obtained by means of Eq. |4]will be de- 
scribed by analytical functions in order to test the fit parame- 
ters for their physical feasibility. An analytical description of 
the spectrum is determined by fitting a series of functions listed 
in Table |2] A ;t^^-min imization alg orithm (utilizing the MINUIT 
package routines, see'James'l998) is employed, starting with the 
first function of the table, a simple power law. The fit is not con- 
sidered valid if the corresponding probability is PfitCf^) < 0.05. 
In this case the next function with more model parameters from 



Table |2] is evaluated. For a given energy spectrum of n data 
points, only functions are examined with « - 1 > degrees of 
freedom. If more than one fit results in an acceptable fit proba- 
bility, an F-Test is used to determine the preferred hypothesis. 
The parameters of the model with more fit parameters are exam- 
ined if the test results in a 95 % probability that description of 
the data has improved. 

3. Exclusion criteria for the EBL shapes 

In the following, arguments to exclude EBL shapes will be pre- 
sented. While the first criteria are based on the expected concav- 
ity of the intrinsic VHE spectra, the second set of criteria relies 
on the integral of the intrinsic VHE emission. 

3.1. Concavity 

Observations have led to the commonly accepted picture that 
particles are accelerated in jets of AGN thereby producing 
non-thermal radiation. The SED of these objects is dominated 
by two components. The first low-energy component from in- 
frared to X-ray energies is due to synchrotron radiation from 
a distribution of relativistic electrons. The second component 
responsible for HE and VHE emission can be explained by 
several different emission models. In leptonic blazar models, 
photons are upscattered by the inverse Compton (IC) pro- 
cess. The involved photon fields originate from synchrotron 
emission (e.g. Bloom & M arscherl 1 19961) . the accretion disk 
('Dermer & Schlick eiser .1993h . or the broad line region (e.g. 
Sikora et al. 1994). In hadronic blazar models, on the other hand, 
y-ray emission is produced either by proton synchrotron radia- 
tion (e.g. Miicke & Protheroe 2001) or photon pion production 
(e.g. iMiicke et al.l i2003). These simple emission models which 
commonly describe the measured data satisfactorily, do not pre- 
dict a spectral hardening in the transition from HE to VHE nor 
within the VHE band. This is also confirmed by observations of 
nearby sources. On the contrary, the spectral slope is thought to 
become softer with energy, either due to Klein-Nishina effects in 
leptonic scenarios and / or a cut off in the spectrum of acceler- 
ated particles. 

However, in more specific scenarios a spectral hardening is 
possible. If mechanisms like, e .g., second order IC sc attering, 
internal photon absorption (e.g. lAharonian et"al]|2008bl) . comp- 
tonization of low frequency radiation by a cold ultra-relativistic 
wind of particles (Aharonian et al. 2002), o r multiple HE and 
VHE y-ray emitting regions in the source dLefaet al.l 12011 ah 
contribute significantly to the overall spectrum, convex curva- 
ture or an exponential pile can indeed occur. Nevertheless, nei- 
ther of these features has been observed with certainty in nearby 
sources. Furthermore, it would imply serious fine tuning if such 
components appeared in all examined sources in the transition 
from the optical thin, i.e. Ty < 1, to optical thick regime, 
Ty > I. This seems unlikely, considering the large number of 
EBL shapes tested. Consequently, EBL shapes leading to an in- 
trinsic VHE spectrum which is not concave will be excluded. 
This expectation is formulated through three test criteria: 

(i) Fermi-LfiJ spectrum as an upper limit. With the launch of 
the Fermi satellite and the current generation of lACTs, there is 
an increasing number of broad-band AGN energy spectra mea- 
sured in the HE and VHE domains. Thus, the least model de- 
pendent approach is to test spectra against a convex curvature in 
the transition from HE to VHE by regarding the spectral index 
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Table 2. Analytical functions fitted to the deabsorbed spectra. 



Description Abbreviation Formula dA'i„t/dis # of parameters 



Power law 

Broken power law 
with transition region 

Broken power law with transition region 
and super-exponential pile up 



PL PL, 2 

BPL PLi X CPLn 4 

SEBPL PLi X CPL,2 X Pile 6 



Double broken power law t^tit^t nr ^nr ^n, er 

^ . . ^ . DBPL PLi X CPL,-, X CPL23 6 
with transition region 

Double broken power law ^^^^^^ PL, x CPL,. x CPL^, x PUe 8 
with super-exponential pile up 11-^. 

Notes. The functions are a power law, a curved power law and a super exponential pile up defined as PL, = Nq E^'^' , CPLjj = 
|l + (ii/iij"^'^"''^^' j , and Pile = exp^(^E / Ep,,^^^^, respectively, where all energies are normalized to 1 TeV. The smoothness parameters / 

are held constant and the break energies E^''"^ are forced to be positive. Only positive pile up, i.e. Ep^^ > 0, is tested. 



measured by Fermi, The, as a limit on the reconstructed intrinsic 
index at VHE, F. Hence, the intrinsic VHE spectrum is regarded 
as unphysical if the following condition is met, 

r + CTstat + CTsys < The - CTHE, stat- (6) 

The statistical error osc^i is estimated from the fit of an analytical 
function to the intrinsic spectrum whereas the systematic uncer- 
tainty (Tsys is used that is estimated by the respective instrumental 
team. The statistical uncertainty ctue. stat of the Fermi-l^PS spec- 
tral index is given by the 2FGL or the corresponding publica- 
tion, see Section m This exclusion criterion will be referred to 
as VHE-HEIndex criterion i n the following. N ote that this is not 
the same criterion as used bv lOrr et"an (1201 ih . They assume that 
the VHE index should be equal to the index measured with the 
Ferm/-LAT. 



DBPL. However, with the VHEConcavity criterion also intrinsic 
spectra can be excluded that show only mild convexity, i.e. no 
exponential pile up. 

3.2. Cascade emission and energy budget 

In this Section two new approaches are introduced that are based 
on the integrated intrinsic emission. These methods rely on a 
number of parameters, whose values are, so far, not accurately 
determined by observations or for which only upper and lower 
limits exist. Therefore, the following two criteria have to be 
regarded as a theoretical motivated possibility to constrain the 
EBL in the future. As it will be shown in Section |5] the final 
upper limits are not improved by these criteria and are, thus, in- 
dependent of the model parameters chosen here. 



(ii) Super exponential pile up. Furthermore, shapes will be ex- 
cluded that lead to an intrinsic VHE spectrum that piles up su- 
per exponentially at highest energies. This is the case if it is 
best described by the analytical functions abbreviated SEBPL or 
SEDBPL, see Table [J] and the pile-up energy is positive within 
a Icr confidence, 

-Epile - CTpiie > 0. (7) 

This additional independent exclusion criterion relies solely on 
VHE observations which are subject to the attenuation in con- 
trast to Fermi-LAiT observations. It will be denoted as PileUp cri- 
terion throughout this study. 

(iii) VHE concavity. In the case that the intrinsic spectrum is 
best described by either a BPL or a DBPL, it is considered as 
convex if the following inequalities are not fulfilled, 

Fi - (Ti < F2 -i- 0-2 
and T2-0-2 < F3 + 0-3 (DBPL), (8) 

and the corresponding EBL shape will be rejected. Again, Icr 
uncertainties of the fitting procedure are used. This criterion will 
be referred to as VHEConcavity. It is very similar to the argu- 
ment formulated in (ii) as intrinsic spectra that show an expo- 
nential rise may often be equally well described by a BPL or 



3.2.1. Cascade emission 

EBL photons that interact with VHE y-rays produce e^e^ pairs. 
These secondary pairs can generate HE radiation by upscatter- 
ing cosmic microwave background (CMB) photons by means 
of the IC process. This initiates an electromagnetic cascade as 
these photons can again und ergo pair production (e.g . ISvenssonI 
1987; Aharonianet al. 1994: iDai et aLll2002t lOoiag et al. 200i 
.Kachelriess et al.ii201 li) . The amount of cascade radiation, that 
points back to the source, depends on the field strength Bigmf 
of the intergalactic magnetic field and its correlation length Ab. 
The values of Bigmf and Ab are unkn own and only upper and 
lower limits exist (see e.g. lNeronov & Semikoz 2009, for a com- 
pilation of limits). If the field strength is large (see Eq. |9]) or if 
the coiTelation length is small compared to the cooling length 
cfcooi of the e^e^ pairs for IC scattering, the pairs are quickly 
isotropized and extended halos of y-ray emission f orm around 
the initial source (e.g Aharonian et al. 1994; Dai et al.l 120021 ; 
lEIviv et ani20 09. Dolag et al. 2009). Furthermore, the time de- 
lay of the cascade emission compared to the primary emission 
depends on Bigmf and Ab. VHE y-ray s need to be produced for 
a sufficiently long period so that the reprocessed radiation is ob- 
servable (e.g. Dermer et al. 201 1). 

The cascade emission has been used to place lower 
limits on Bigmf and Ab by assuming a certain EBL 
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model (iNeronov & Vovkl |20 1 Ot iTavecchio etani2010ll2011 



iDermer et alJ 120111: iDolagetalJ 120111: lHuanetairi201lh 
Conversely, one can place upper limits on the EBL density under 
the assumption of a certain magnetic field strength. This novel 
approach is followed here whereas, in previous studies, the cas- 
cade emission is neglected when deriving upper limits on the 
EBL density. A higher EBL density leads to a higher produc- 
tion of e^e^ pairs and thus to a higher cascade emission that is 
potentially detectable with the Fen«/-LAT. If the predicted cas- 
cade radiation exceeds the observations of the Fermi-LAI, the 
corresponding EBL shape can be excluded. Conservative upper 
limits are derived if the following assumptions are made: (i) The 
HE emission of the source is entirely due to the cascade, (ii) 
The observed VHE spectrum is fitted with a power law with a 
super exponential cut off at the highest measured energy of the 
spectrum. This minimizes the reprocessed emission and allows 
to consider only the first generation of the cascade, (iii) The e^e^ 
pairs are isotropized in the intergalactic magnetic field, minimiz- 
ing the reprocessed emission. This condition is equal to the de- 
mand that the deflection angle ■& of the particles in the magnetic 
field is a! t:. Assuming Ab » cfcooi, the deflection angle for elec- 
trons with an energy ymc^ » E/2, where £ is the energy of the 
primary y ray, can be approximated by (ITavecchio et al.ll2010l: 
INeronov & Vovgi20Toh 



Rl 1 10-15 g/ \1Q(>) 



(9) 



with Zr the redshift where the IC scattering occurs and /?l 
the Larmor radius. The IC scattered e^e" pairs give rise to y- 
rays with energy e ~ j^Iivcmb ~ 0.63(£'/TeV)^ GeV, with 
hvcMB - 634 ^leV the peak energy of the CMB. The y fac- 
tor in Eq. |9| can be eliminated in favor of e, and, solving for 
fiiGMF, the pairs are isotropized if Bigmf ~ 4.2 x 10"'^ (1 + 
Zr)\e/GeY) G ^ 5 x 10"'^ G for e = 100 GeV, the maximum 
energy measured with the Fermi-LAT considered here and the 
maximum redshift where the IC scattering can occur, i.e. the red- 
shift of the sourceQ This value of Z^igm f is in acco r dance with 
all experimental bounds (see e.g. iNeronov & Semikozl l2009l 
especially Figures 1 and 2) For correlation lengths Ag » 
cfcooi ~ 0.65(£/TeV)-i(l + z,)-'*Mpc ^ (9(Mpc) the most 
stringent constraints come from Faraday rotation measurements 
(Kronberg & Simard-Normandin 1976; Blasi et al. 1999) which 
limit BiGMF ^ 10"^ G. Furthermore, the adopted value cannot 
be excluded neither with possib le observations o f deflections of 
ultra-high energy cosmic rays (iLee etalJll99l nor with con- 
strained simulations of magnetic fields in gal axy clusters, both 
setting an upp er limit on Bigmf ^ 10 '^G ( Dolag et al.ll2005t 
iDonnert et all2009l) . 

For value of Bigmf, the cascade emission is detectable if a 
steady y-ray emissi on of the source for the last At ^ 10^ year s 
is assumed (see e.g. lDermer et alj[201 ll : iTkvecchio et al.ll201 ih . 
Other energy loss channels apart from I C scattering like syn- 
chrotron radiation or plasma instabilities dBroderick et al.ll201 ll) 
are neglected. However, if the latter are present, the field strength 
is even higher, or the lifetime of the VHE source is shorter, no 
significant cascade emission is produced or it has not reached 
earth so far 

The cascade emis sion F(e) is calculated wit h Eg. IB .3 1 in 
Appen dixlBlfollowing lTavecchio et sl\ (1201 ih and'Perm er et alJ 
(1201 lb . For isotropy, the observed cascade emission has to be 
further modified with the solid angle Q.^ ~ i'^'^o which the 



' Accordingly, this S-field value ensures isotropy regardless were 
the IC scattering occurs. 



intrinsic blazar emission is collimated where 0c is the semi- 
aperture of the irradiated cone. For blazars one has Gc ~ I /Fl 
where Fl is the bulk Lorentz factor of the plasma of the jet. T he 
observed emission is then found to be (ITavecchio et al.ll201 ll) 



(10) 



where the factor of two accounts for the contribution of both jets 
in the isotropic case. The exclusion criterion for an EBL shape 
at the 2cr level reads 



£obs(^meas) ^ ^meas "t" 2crjy 



(11) 



where e^eas, ^^meas, o",j,ea,s are the measured energy, flux and sta- 
tistical uncertainty reported in the 2FGL, respectively. In the 



case that the source is not detected F„ 



= and cr„ 



rep- 



resents the 1 cr upper limit on the flux. As an example. Figure 
|2| shows the observed and intrinsic VHE spectrum for a spe- 
cific EBL shape of flie blazar lES0229-t- 200 together with the 
Fermi upper limits (ITavecchio et alJl2010l) . The different model 
curves demonstrate the degeneracy between the different param- 
eters entering the calculation. The EBL shape used to calculate 
the intrinsic VHE spectrum is not excluded in the isotropic case 
since the emission does not overproduce the Fermi upper limits. 
This is contrary to the case of Bigmf - 10^^^ G and Af = 3 years 
where the predicted cascade flux exceeds the Ferm/-LAT upper 
limits in the 1 - 1 GeV range. To obtain conservative upper limits 
of the EBL, only the isotropic case is assumed in the following, 
i.e. BiGMF = 5 X lO"'-' G which implies that the source has to be 
steady for a lifetime of At > 10^ years. Furthermore, a Lorentz 
factor of Fl - 10 is generically assumed for all sources. 

3.2.2. Total energy budget 

The jets of AGN, the production sites for HE/VHE emission, 
are powered by the accretion of matter onto a central black hole 
(e.g. Urrv & Padovani 1995). If the radiation escapes isotropi- 
cally from the black hole, the balancing of the gravitational and 
radiation force leads to the maximum possi ble luminosity due 
to acc retion, the Eddington luminosity, (e.g. iDermer & MenonI 
I2009h 

iedd(M.) a 1.26 X 10^** ^ ergs s"' 



Mr, 



(12) 



where M, is the black-hole mass normalized to the mass of 
the sun, Mq. Assuming that the total emission of an AGN is 
not super-Eddington, the Eddington luminosity is the maximum 
power available for the two jets, Pjet < Ledd/2, which is a sum 
of several contributions which all can be represented as (e.g. 
ICelotti & Fabiarilll993l: iBonnoh et alJl201 Ih f, = ttR'^-F^/ScU'i, 
in the case that the radiation is emitted by an isotropically radi- 
ating relativistic plasma blob in the comoving frame. The blob 
of radius R' in the comoving frame moves with a bulk Lorentz 
factor Fl and corresponding speed /3c, U'. is the comoving en- 
ergy density. The energy density of the produced radiation is 
t/; = L'liAnR'^c) = LliAnS'l^R'^c). The last equality con- 
nects the comoving luminosity with the luminosity in the lab 
frame via the Doppler factor 5d = [(1 + z)TiX^ -jScos^)]"' ^ 
2Fl/[(1 + z)(l + 6^T\)] where 6 is the angle between the jet 
axis and the line of sight. The approximation holds for <K 1 
and Fl » 1. Assuming 6 x 6^, the Doppler factor and the bulk 
Lorentz factor are equal up to the redshift factor \ + z. So ~ T]^. 
The power produced in rad iation is a robust lo wer limit for the 
entire power of the jet (e.g. lBonnoli et al.ll201 lb . 



A- ^ LliATl) < Pjet < Ledd/2 



(13) 
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Fig. 2. Cascade emission for a certain EBL shape and the VHE 
specti-um of lES 0229+200 (Ahai-onian et al. 2007a). The ob- 
served spectrum (dark red points and Une) is fitted with a power 
law with an exponential cut off and corrected for the EBL ab- 
sorption (dark blue dashed line and points). The green lines show 
the cascade emission resulting from the reprocessed flux (light 
gray shaded area) for a constant emission over the last three 
years and different magnetic field strengths. The red dotted line 
shows the reprocessed emission if the e^e^ pairs are isotropized. 
The latter does not overproduce the Fermi upper limits (black 
diamonds iTavecchio et al.l |2010) and hence the corresponding 
EBL shape is not excluded. The light and dark gray area together 
are equal to the integrated flux that is compared to the Eddington 
luminosity (see Section [3.2.2b . 



Solving the inequality for the observed luminosity, one arrives 
at an additional exclusion criterion for EBL shapes, namely, if 
the intrinsic energy flux at VHE is larger than the associated 
Eddington energy flux. 



(1+z) 



2-T; 



d£ 



-dE > 



(14) 



where and E^^y, are the minimum and maximum energy 
of the intrinsic VHE spectrum which is described with a power 
law with index Pint. The the factor (1 + zf'^^'"' accounts for the 
K-correction and c/l is the luminosity distance given by 



ii +z)c r 

Ho j 



dz' 



(15) 



For a conservative estimate, M, + ctm. is used in the calcula- 
tion of Ledd- The assumption of a non super-Eddington lumi- 
nosity is, however, somewhat speculative as super-Eddington 
emission has been observed e.g. in the variable source 3C 454.3 
(Abdo et al. 2011b). In Section |5] it will be shown that the ca- 
pability of the Eddington criterion to exclude EBL shapes is ex- 
tremely limited. For this reason, only steady sources (listed in 
TableUJi will be considered for this criterion. Here, it is only em- 
phasized that it is in principle possible to constrain the EBL with 
this argument. 

Excluding EBL shapes with cascade emission (Eq. fTTT i and 
the total energy budget of the source (Eq. [T4l i. will be referred 



to as the IntVHELumi (short for intrinsic VHE luminosity) crite- 
rion. 

4. VHE AGN Sample 

In the past four years, the number of discovered VHE emit- 
ting AGN has doubled. In this section samples of VHE spec- 
tra are defined that are evaluated with the VHE-HEIndex, 
PileUp and VHEConcavity criteria (Section l4.ll) and with the 
IntVHELumi criterion (Section l4~2b . 

4. 1 . Sample tested against concavity criteria 

For this part of the analysis, 22 VHE spectra from 19 differ- 
ent sources are used. AGN are included in the sample only if 
their redshift is known, there is no confusion with other sources, 
and they are detected with the Fen7i/-LAT. This excludes 
the known VHE sources 3C66A and 3C66B, lES 0229+200, 
PG 1553 + 113, and S50716+714. Two spectra from the same 
source are only considered if they cover different energy ranges. 
Furthermore, the radio galaxies Centaurus A and M 87 are not 
included since they are too close and measured at energies that 
are too low to yield any constraints of the EBL density. Spectra 
that are a combination of several instruments are not included 
due to possible systematic uncertainties. If two or more spectra 
are available for a variable source, the VHE spectrum is cho- 
sen that is measured simultaneously with Ferm/-LAT observa- 
tions. If the Fermi spectrum is best described with a logarithmic 
parabola, the spectral index determined at the pivot energy is 
used for the comparison with the intrinsic VHE spectra. The en- 
tire AGN sample is listed in Table |3] together with the redshift, 
the energy range, the spectral index at VHE energies, the index 
measured with the ferm/-LAT, the variability index given in the 
2FGL, and the corresponding references. 

AGN are known to be variable sources both in overall 
flux and spectral index. This poses a problem for the VHE- 
HEIndex criterion as it relies on the comparison of Fermi- 
LAT and lACT spectra. To address this issue, one can roughly 
divide the overall source sample into three categories: 

1 . Steady sources in the Fermi-LAr energy band. In this cat- 
egory, all sources are assembled that show a variability in- 
dex < 41.64 in the 2FGL which corresponds to a likelihood 
fit probability of more than 1 % th at the source is steady 
dThe Fermi-LAT collaboration! 1201 ll) . For these sources si- 
multaneous measurements are not required regardless if they 
are steady (like lES 1101-232, Aharonian et al. 2007b) or 
variable (like H 1426+428, see below) at VHE. This does not 
affect the upper limits derived here because the Fermi index 
remains valid as a lower limit independent of the VHE index. 
These sources are marked as "steady" in the last column of 
Tables 

2. Variable sources with simultaneous measurements. Some 
of the variable sources were observed simultaneously 
with the Ferm;-LAT and lACTs in multi wavelength cam- 
paigns, namely PKS 2155-304 with HESS (Aharonian et alj 
2009), and PKS 1222+21 with MAGIC ( Aleksic et al. 201 C 
Instead of the spectral slopes given in the 2FGL, the Fermi- 
LAT spectra from these particular observations are used to 
test the EBL shapes. These sources are marked as "simul" in 
the last column in Table |3] Note, however, that the observa- 
tion times might not be equal for the individual instruments 
since the sensitivities for the Fermi-LAI and lACTs are dif- 
ferent. Nevertheless, the arising systematic uncertainty is 
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Table 3. VHE AGN spectra used in this study. If not stated otherwise in the text, the Fermi slope and variability index are taken 
from the 2FGL. 



Source 


Redshift 


Experiment 


Energy Range 
(TeV) 


VHE Slope 

r ± 0"siai ± 0".sys 


Fermi Slope 
r ± o-su,, 


Variability Index 


Reference 


Comments 


Mkn421 


0.031 


HESS 


1.75 - 23.1 


2.05 ± 0.22 


1.77 


± 


0.01 


112.8 


(1) 


hardest index 


Mkn501 


0.034 


MAGIC 


0.17-4.43 


2.79 ± 0.12 


1.64 


± 


0.09 


72.33 


(2) 


hardest index 


Mkn501 


0.034 


HEGRA 


0.56-21.45 


1.92 ± 0.03 ± 0.20 


1.64 


± 


0.09 


72.33 


(3) 


hardest index 


lES 2344+5 14 


0.044 


MAGIC 


0.19-4.00 


2.95 ± 0.12 ± 0.20 


1.72 


± 


0.08 


28.13 


(4) 


steady 


Mkn 180 


0.045 


MAGIC 


0.18- 1.31 


3.25 ± 0.66 


1.74 


± 


0.08 


19.67 


(5) 


steady 


lES 1959+650 


0.048 


HEGRA 


1.52- 10.94 


2.83 ± 0.14 ± 0.08 


1.94 


± 


0.03 


52.30 


(6) 


hardest index 


lES 1959+650 


0.048 


MAGIC 


0.19-2.40 


2.58 ± 0.18 


1.94 


± 


0.03 


52.30 


(7) 


hardest index 


BIj Lacertae 


0.069 


MAGIC 


0.16-0.70 


3.6 ± 0.5 


2 11 






267.0 


(8) 


hardest index 


PKS 2005-489 


0.071 


HESS 


0.34-4.57 


3.20 ± 0.16 ± 0.10 


1.90 


+ 


0.06 


68.86 


(9) 


hardest index 


RGB J0152+017 


0.080 


HESS 


0.31-2.95 


2.95 ± 0.36 ± 0.20 


1.79 


± 


0.14 


27.73 


(10) 


steady 


PKS 2155-304 


0.116 


HESS 


0.25-3.20 


3.34 ± 0.05 ± 0.1 


1.81 


± 


0.11 


262.9 


(11) 


simul 


RGB J07 10+591 


0.125 


VERITAS 


0.42-3.65 


2.69 ± 0.26 ± 0.20 


1.53 


± 


0.12 


29.86 


(12) 


steady 


H 1426+428 


0.129 


HEGRA 


0.78-5.37 




1.32 


± 


0.12 


22.16 


(13) 


steady 


IBS 0806+524 


0.138 


MAGIC 


0.31-0.63 


3.6 ± 1.0 ± 0.3 


1.94 


± 


0.06 


37.80 


(14) 


steady 


H 2356-309 


0.165 


HESS 


0.23-1.71 


3.06 ± 0.15 ± 0.10 


1.89 


± 


0.17 


20.19 


(15) 


steady 


lES 1218+304 


0.182 


MAGIC 


0.09 - 0.63 


3.0 ± 0.4 


1.71 


± 


0.07 


40.00 


(16) 


steady 


lES 1218+304 


0.182 


VERITAS 


0.19-1.48 


3.08 ± 0.34 ± 0.2 


1.71 


± 


0.07 


40.00 


(17) 


steady 


IBS 1101-232 


0.186 


HESS 


0.18-2.92 


2.88 ± 0.17 


1.80 


± 


0.21 


25.74 


(18) 


steady 


lES 1011+496 


0.212 


MAGIC 


0.15-0.59 


4.0 ± 0.5 


1.72 


± 


0.04 


48.05 


(19) 


hardest index 


lES 0414+009 


0.287 


HESS 


0.17-1.13 


3.44 ± 0.27 ± 0.2 


1.98 


± 


0.16 


15.56 


(20) 


steady 


PKS 1222+21 " 


0.432 


MAGIC 


0.08-0.35 


3.75 ± 0.27 ± 0.2 


1.95 


± 


0.21 


13030 


(21) 


simul 


3C279 


0.536 


MAGIC 


0.08 - 0.48 


4.1 ± 0.7 ± 0.2 


2.22 


± 


0.02 


2935 


(22) 


hardest index 


lES 0229+200" 


0.140 


HESS 


0.60-11.45 


2.5 ± 0.19 ± 0.10 










(23) 





Notes. See the text for details on the Comments column. There was no simultaneous measurement during the 0.5 h in which MAGI C detected the 
sourc e. However, the index used in high energies was extracted from Fermi data in the 2.5 h before and after the MAGIC observation lAleksic et all 
1201 Ih . The spectrum is only tested against the IntVHELumi criterion as it is not detected with the Ferrai-LAT. 



Referen ces. (1) [T luczykonj |2011h; (2) [Abdo et al i (1201 lal); (3) lAharo nian et al. (1999); (4) Albert et al ( 2007a): (5) [Albert et al. I 
(120069); (6) lAharonian et all J2003bh; (7) Tadiaferri et al.l (l2008l); (8) Albert et al. (2007b ); (9) ,HESS Colla boration et all (l2010b(); 
(10) Aharonian et al. (20 0ji); (1 1) lAharonian e t al. (2009); (12) Acciari et al. (2010); (13) Aharoni aiTet al.l (l2003al); ( 14) Acciari et al. (2009^ 
(15) ,HESS Col laboration^al.' ("20 1 Oah; (16)"lAlbert et al . (2006a); (1 7) .Acciari et al. (2009b); (18) lAharonian et al.l 1 2006); (19) .Albert et 5] 
( I2007d) ; (20)lThe HESS Collabor ation et alj ( 120121) ; (21)[Aleksic e t al] ( 120111) ; (22) [Albert et al] ilOO^ : (23) lAharonian eTaLr ( l2007al) 



negligible for the sources under consideration. In the case of 
PKS 2155-304, the source was observed in a quiescent state 
where no fast flux variability is expected. PKS 1222-1-21, 
on the other hand, was observed in a HE flaring state and 
Ferwi-LAT observations are not avail able the 30 minutes of 
MAGIC observations. Instead, Ale ksic et al .l (|2011) derive 
the Fermi spectrum from 2.5 hrs of data encompassing these 
30 minutes. This is justified, since the source remained in 
this high flux state for s everal days with little spectral vari- 
ations (cf. Figure 2 in Tanakaet al. 2011). Accordingly, the 
maximum time lag allowed for observations to be considered 
as simultaneous is of the order of an hour. 
3. Variable sources not simultaneously measured. For some 
variable sources, no simultaneous data are available, namely, 
lES 1011+496, lES 1959+650, 3C279, BLLacertae, the 
flare spectra of Mkn 501 and Mkn 421, and PKS 2005-489 
(see Table [3] for the references). In these cases, the litera- 
ture was examined for dedicated Fermi-LAiT analyses of the 
corresponding sources in order to find hardest spectral index 
pubHshed. In the cases of lES 1011+496, lES 1959+650, 
PKS 2005-489 and Mkn 421 the indices reported in the 
2FGL are the hardest published so far. The har dest indices 
for BL Lacertae and Mkn 501 are obtained by Abdo et al.' 
(l20I0al) and Abdo et al. (201 la), respectively, see Table[3]for 
the corresponding values. The distant quasar 3C 279 was ob- 
served with the Fermi-LAL? during a y-ray flare in 2009 and 
the measured spect ral indices vary be tween ~ 2 and ~ 2.5 
(compare Fig. 1 in Abdo et al. 2010c'). Thus, the catalog in- 
dex of 2.22 + 0.02 is appropriate to use. Table prefers to 



all the spectra discussed here as "hardest index" in the last 
column. 

Additional uncertainties are introduced for the VHE observa- 
tion with a maximum time lag between the measurement and the 
launch of the Fermi satellite, which is the case for Mkn 501 and 
H 1426+428. In the case of H 1426+428, no detection has been 
reported after the HEGRA measurement in 2002 at VHE which 
might suggest that the source is now in a quiescent state. The 
2002 spectrum with an observed spectral index of F = 1 .93+0.47 
is used in this study. The hard spectrum promises stronger limits 
with the VHE-HEIndex criterion than the 2000 spectrum which 
has a spectral slope of F = 2.79 ± 0.33. The source showed 
a change in flux by a factor of 2.5 between the 1999/2000 
and 2002 observation runs but the spectral slope remained con- 
stant (Aharonian et al. 2003a). Additionally, the Fermi index of 
1 .32 + 0.12 is the hardest of the entire sample and, in summary, 
it is chosen to include the source in the study. As for Mkn 501, 
the spectrum of the major outbreak was measured up to 21 TeV 
and, consequently, it is a promising VHE spectrum to constrain 
the EEL density at FIR wavelengths. As it turns out, it excludes 
most shapes due to the PileUp and VHEConcavity criteria. These 
criteria are independent of the Fermi index and, therefore, not af- 
fected by the diff'erence in observation time. 

4.2. Sample tested against intrinsic VHE luminosity 

For the integral criterion presented in Section 13.2.11 and 13.2.21 
only spectra from steady sources are used in order to avoid 
systematic uncertainties introduced by variability. Only spec- 
tra are examined which suffer from large attenuations and are 
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Table 4. Sources used to exclude EBL shapes with the 
IntVHELumi criterion. 



Black-hole mass 



lES 0229+200 


9.16 


+ 


0.11 


lES 0414+009 




9.3 




lES 1101-232 




9 




lES 1218+304 


8.04 


+ 


0.24 


H 1426+428 


8.65 


+ 


0.13 


H 2356-309 


8.08 


+ 


0.23 


RGBJ0152+017 




9 




RGBJ0710+591 


8.25 


+ 


0.22 



Notes. The black hole masses M. are taken from IWaene j i2008h ex- 
cept from RGBJ0710+591 a nd I BS 0414+00 9 for which the masses 
are given in I Woo et al.l ( l2005h and lUrrv et al.l (2000), respectively. No 
measurements of the central black hole masses of lES 1101-232 and 
RGB J0152+017 are available so the fiducial value of M. = 10' Mq is 
used here. 

measured at energies beyond several TeV. These spectra are 
the most promising candidates for constraints as they show the 
highest values of integrated intrinsic emission. On the other 
hand, the spectrum of 1E S 022 9+200 which has been rean- 
alyzed by iTavecchio et alj (1201 0) can be tested against the 
IntVHELumi condition as upper limits on the HE flux suffice and 
no spectral information is required for this criterion. Otherwise, 
the same selection criteria apply as for the sample tested against 
concavity (known redshift, etc.). The VHE spectra evaluated 
with the IntVHELumi criterion together with the central black- 
hole masses of the corresponding sources are summarized in 
Table m 

5. Results 

The upper limits on the EBL density are derived by calculating 
the envelope shape of all allowed EBL shapes. The influence 
of the different exclusion criteria is examined by inspecting the 
envelope shape due to the VHE-HEIndex argument alone and 
successively adding the other criteria and reevaluating the re- 
sulting upper limits. Furthermore, the impact of the VHE spec- 
tra responsible for the most stringent limits in the optical, MIR 
and FIR will be investigated by excluding these spectra from the 
sample and inspect the change in the upper limits. 

Figure [3] shows a histogram of the fractions of rejected 
shapes by each VHE spectrum, where the different colors rep- 
resent the different criteria that lead to the exclusion of an EBL 
shape. It should be noted that individual shapes can be re- 
jected by several criteria at the same time, and, therefore, the 
different columns may add up to a number larger than indi- 
cated by the total column. Results for spectra that exclude no 
(BLLacertae, lES 2344+514, and MknlSO) or less than 10% 
of an EBL shapes (the MAGIC spectrum of lES 1959+650, 
the HESS spectra of RGB J0152+017 and Mkn 421, as well as 
the HEGRA spectrum of lES 1959+650) are not shown. Most 
EBL shapes are excluded by the VHE spectra of H 1426+428, 
lES 1101-232, and Mkn 501. The influence of H 1426+428, and 
Mkn 501 on the limits in the MIR and FIR and of 3C279 to- 
gether with PKS 1222+21 in the optical will be examined by 
excluding these spectra from the sample. These sources provide 
strong constraints in the respective wavelength bands. Note that 
removing lES 1 101-232 from the source sample does not change 



the upper limits since a number of spectra of sources with com- 
parable redshifts (cf. Table |3]l exclude the same EBL shapes as 
lES 1101-232, e.g. lES 0414+009, the VERITAS spectrum of 
lES 1218+304, H 2356-309, and PKS 2005-489. 

Different combinations of exclusion criteria are shown in the 
panels of Figure |4] Each panel depicts the limits for the com- 
plete spectrum sample and, additionally, the resulting EBL con- 
straints if the spectra discussed above are omitted. By itself, the 
VHE-HEIndex criterion gives strong upper limits in the optical 
and MIR on the EBL density if all spectra are included (upper 
left panel of Figure HJ. In the optical, the limits are dominated 
by the spectra of 3C279 and PKS 1222+21, so, consequently, 
the restrictions are significantly weaker without these spectra 
(dashed line in Figure |4]i. The spectra are influenced most by 
changes of the EBL density in the optical which is inferred from 
the maximum energies of 480 GeV and 350 Ge V for 3C 279 and 
PKS 1222+21, respectively. They translate into maximum cross 
sections for pair production at 0.6 ^im (3C279) and 0.43 ^im 
(PKS 1222+21), see Eq.|3] The constraints are almost unaltered 
if only one of these spectra is excluded from the sample. In the 
MIR, the spectrum of H 1426+428 provides firm limits on the 
EBL density whereas scarcely any EBL shape is rejected due to 
the spectrum of Mkn 501 with the VHE-HEIndex criterion. 

The combination of the VHE-HEIndex and VHE- 
Concavity criterion strengthens the upper limits between 
2 [,im and 10 [^im, as shown in the upper right panel of Figure 
|4] Convex intrinsic spectra are the result of an EBL density 
with a positive gradient between lower and higher wavelengths 
and, thus, a combination with the VHE-HEIndex criterion is 
necessary to exclude shapes with a high EBL density that are 
rather constant in wavelength. Therefore, on their own, neither 
the PileUp nor the VHEConcavity criterion provide strong upper 
limits. Combining the PileUp and VHE-HEIndex arguments 
results in very similar limits as the combination of the VHE- 
HEIndex and VHEConcavity criterion. This degeneracy between 
the PileUp and VHEConcavity criterion is also demonstrated 
in Figure |5] The spectrum of Mkn 501 corrected with a certain 
EBL shape shows a strong exponential rise at highest energies 
but is best described with a double broken power law. The 
combination of the PileUp and VHEConcavity together with the 
VHE-HEIndex criterion yields robust upper limits in the FIR as 
displayed in the lower right panel of Figure |4] The constraints 
in the FIR are entirely due to the spectrum of Mkn 501 although 
the spectrum of Mkn 421 is also measured beyond 20 TeV and 
both sources have a comparable redshift. However, the spectrum 
of Mkn 421 rejects far less shapes than Mkn 501. Indeed, an 
exponential rise is observed in intrinsic spectra of Mkn 421 for 
certain EBL realizations (e.g. the corrected Mkn 421 spectrum 
in Figure |5j but a power law is found to be the best description 
of the spectrum. 

Compared to the VHE-HEIndex criterion alone, the com- 
bination with the IntVHELumi criterion leads to improved up- 
per limits only if H 1426+428 is discarded from the sample 
(lower left panel of Figure |4|. Most shapes are rejected by the 
VHE spectrum of lES 0229+200 which is also the sole spec- 
trum which excludes a very limited number of shapes with the 
Eddington luminosity argument. Remarkably, the spectrum ex- 
cludes more than 60 % of all shapes. The IntVHELumi criterion 
has the most substantial effect in the infrared part of the EBL 
density as the highest energies of the spectrum of lES 0229+200 
contribute most to the integral flux. The maximum energy mea- 
sured in the spectrum is 1 1 .45 TeV and thus the limits are most 
sensitive to changes in the EBL around 14 (im. The influence 
of the choice of the bulk Lorentz factor Fl (and hence of the 
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Fig. 3. Histogram of the fraction of excluded shapes of the different VHE spectra. The columns show the total fraction of rejected 
shapes as well as the fraction excluded by the different criteria. The column labeled "Curvature" combines the VHEConcavity and 
PileUp criteria. Spectra that allow more than 90 % of all shapes are not shown. 



Doppler factor 5d since Fl ~ 5d is assumed) on the envelope 
shape can be seen from Figure |6l where the upper limits are 
shown for Fl = 5, 10, and 50. As Fl enters quadratically into 
the calculation of the flux (cf. Eq. [TOb and for the Eddington lu- 
minosity (Eq. [T4b the choice of the value of Fl is critical for 
the number of rejected EBL shapes. The bulk Lorentz factor is 
unknown for the sources tested with IntVHELumi and for the 
combination with the other criteria Fl - 10 is generically cho- 
sen. However, even with this oversimplified choice of Fl, the 
IntVHELumi criterion does not lead to improvements of the up- 
per limits compared to the combination of the VHE-HEIndex , 
VHEConcavity and PileUp criteria. Conversely, this implies 
that the final upper limits will not depend on the specific choice 
of model parameters and assumptions that enter the evaluation 
of the IntVHELumi criterion. 

The final result for the upper limits is the combination of 
all criteria and all VHE spectra, shown in Figure [T] It is the 
envelope shape of all allowed EBL realizations, cf. Figure |7^, 
which itself is excluded by several VHE spectra and it should, 
thus, not be regarded as a possible level of the EBL density. 
For the maximum energy of all VHE spectra of 23.1 TeV, the 
cross section for pair production peaks at a wavelength of the 
EBL photons of /l» ^ 29 More than half of the interactions 
occur in a narr ow interval A/1 - (1 + 1/2) A« around the peak 
wavelength (e.g'Aharo nian et alJ2006l) and hence the constraints 
are not extended beyond 100 [im. Albeit including the evolution 
of the EBL with redshift, the derived upper limits are below 
5 nWm"^ sr"' in the range from 8 (xm to 31 (xm. A comparison 
of the constraints with previous works is shown in Figure |7j). 
Above 30 ^im, the constraints are consistent with those derived 
in MR07. For wavelengths betwe en 1 \m\ and 4 >xm the li mits are 
in a ccordance w ith the results of lAharonia n et al. ( 2006, 2007a') 
and lAlbert et"al] (l2008i) . The strong hmits of Albert et al.. (2008^ 
who utiUzed the spectrum of 3C 279 are not reproduced. Note, 
however, that these limits are derived by changing certain free 
parameters (e.g. fraction of UV ernission escaping the galaxies) 
of the EBL model of iKneiske et al.l (120021) while the current ap- 
proach allows for generic EBL shapes. Consequently, an EBL 
shape with a high density at UV / optical wavelengths followed 
by a steep decline towards optical and NIR wavelengths pro- 
duces a soft intrinsic spectrum of 3C 279 that cannot be excluded 
by any criterion. Furthermore, an inspection of the spectrum of 
3C 279 shows that the fit will be dominated by the first two en- 



ergy bins due to the smaller error bars. Thus, a convex spectrum 
is often still sufficiently described with a soft power law. In gen- 
eral, it should be underlined that all of the above limits from re- 
cent studies use a theoretically motivated bound on the intrinsic 
spectral slope of F = 1.5. 

The upper limits derived here are not in confl i ct wit h 
the EBL mode l calculations of 'Franceschini et al.) ( l2008l) . 
i Kneiske & Dold (|2010) and Dommguez et al. (201ib and are 
compatible with the lower limi t galaxy nu mber counts de- 
rived from Spitzer measurements dFa zio et al]|2004i) (see F igure 
17}:). In the FIR, t he mo dels of iFranceschini et all (l2008l) and 

Ido 



iDomfnguez et al.1 (12011 ') lie above the derived upper limits, 
though one should note that the EBL limit at these wavelengths 
relies on a single spectrum (Mkn501). Between ~ 1 (im and 
~ 14 ^.m there is, however, a convergence between the upper 
limits and model calculations and at 13.4(xm the EBL is con- 
strained below 2.7 nWm"^ sr"', just above the EBL models. 
This leaves not much room for additional components such as 
Populati on III stars and implies that the direct measurements of 
Matsumoto et alJ ( 20051) are foreground dominated as discussed 
m iDwek et alJ ( l2005h . The EBL models, the upper limits from 
previous works, and the results derived here are shown together 
in Figure [8] 

In general, most of the tested EBL shapes are excluded by 
more than one spectrum (Figure |9]l. While 0.23 % of all EBL 
shapes are excluded by only one of the spectra in the sample, 
the majority of shapes (93 %) is rejected by five spectra or more. 
Figure |7}l shows the limits for different minimum numbers of 
VHE spectra that rule out an EBL shape. From NIR to MIR 
wavelengths, the limits are only slightly worsened if at least two 
spectra are required to exclude EBL shapes. If at least five spec- 
tra are ought to reject an EBL shape, the EBL density remains 
confined below 40 nW sr ' m"^ in the optical. Thus, from opti- 
cal to MIR wavelengths, the hmits are robust against individual 
spectra that possibly have a peculiar intrinsic shape due to one of 
the mechanisms discussed in Section lTTI Especially in the MIR 
and FIR, however, the limits are weakened as they mainly de- 
pend on two spectra, H 1426-1-428 and Mkn501. This underlines 
the need for more spectra measured beyond several TeV in order 
to draw conclusions about the EBL density in the MIR and FIR 
from VHE blazar measurements. 

Given the similarities in procedures used, the systematic un- 
certainties of the limits derived in this study are similar to the 
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Fig. 4. Limits on the EBL density for different exclusion criteria. The solid line is the envelope shape of all allowed shapes from 
the combination of all VHE spectra whereas the dashed curve shows the envelope shape without considering the VHE spectrum 
of 3C279. The dashed dotted line is the envelope shape without H 1426+428 and the dotted Une displays the upper limits without 
MknSOl. 



ones derived in MR07 and . They have been estimated to be 3 1 % 
in optical to near infrared and 32 - 55% in mid to far infrared 
wavelengths, mainly from the grid spacing and the uncertainties 
on the absolute energy scale of gr ound based VHE instruments 
which is taken to be 15 %. Note that lMever et alj (1201 0) achieved 
a cross-calibration using the broadband SED of the Crab Nebula 
between the Ferm/-LAT and lACTs by shifting the lACT mea- 
surements by ~ 5 % in energy. As it turns out, additional uncer- 
tainties arise from the phenomenological description of the EBL 
evolution (< 4 % for a redshift z - 0.2 and < 10 % for z = 0.5, 
iRaue & Mazinl2008l) . Uncertainties in the calculation of the cas- 
cade emission are caused by the choice of the model parameters 
which are, however, difficult to quantify. The same applies for 
the assumption that steady sources do not show super-Eddington 
luminosities. Since the most stringent limits do not rely on these 
exclusion criteria, these uncertainties do not affect the final re- 
sults of the upper limits. 

Additionally, the measurement capabilities of the Fermi- 
LAT affect the VHE-HEIndex criterion and hence the upper lim- 
its. While the 2FGL does not quote the systematic errors on the 
individual spectral indices, it gives a number of sources of sys- 
tematic errors: the eff'ective area, the diff'use emission model, and 
the handling of front and back converted events. The systematic 
error on the eff'ective area is estimated to be between 5 % and 



10 %, while the errors on the diff'use emission model mainly ef- 
fects sources inside the galactic plane. Furthermore, the isotropic 
emission for front and back converted events is assumed to be 
equal. This leads to underestimation of the flux below 400 Me V 
and might produce harder source spectra. As harder spectra in 
the Fermi-hPil band weaken the upper limits, the results derived 
here can, again, be regarded as conservative. 



6. Summary and Conclusions 

In this paper, new upper limits on EBL density over a wide wave- 
length range from the optical to the far infrared are derived, uti- 
lizing the EBL attenuation of HE and VHE y-rays from distant 
AGN. A large number of possible EBL realization is investi- 
gated, allowing for possible features from, e.g., the first stars. 
Evolution of the EBL density with redshift is taken into account 
in the calculations using a phenomenological prescription (see 
e.g. Rau e & Mazinll2008h . A large sample of VHE spectra con- 
sisting of 23 spectra from 20 different sources with redshifts 
ranging from z - 0.031 to 0.536 is used in the analysis. The 
VHE spectra are corrected for absorption and subsequently in- 
vestigated for their physical feasibility. Two basic criteria are 
examined: (1) concavity of the high energy part of the spec- 
trum spanning from HE to VHE and (2) total integral flux in 
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Fig. 7. Upper limits derived in this study, (a) The envelope shape (upper limits) of all allowed shapes (dark gray Unes). Also 
shown are the grid points as lig ht gray bullets, (b) The constraints compared to the the upper limits of MR07, Aharonian et aJJ 
()2006, 2007a), and Albert et aj j (, , 2008). (c) U pper limits of this study together with three EBL models (Franceschini et al. 2003^ 
[Kneiske & D o"Sl2010r iDomfnguez et alJl201 ih . (d) Upper limits requiring different minimum numbers of VHE spectra that exclude 
an EBL shape. 



the VHE, a novel way to probe the EBL density. For the for- 
mer criterion, spectra from the ferm/-LAT at HE are used as 
a conservative upper limit, combined with criteria on the over- 
all VHE concavity. This is a more conservative argument than 
a theoretically motivated bound on the intrinsic spectral index 
at VHE of, say, F = 1.5. This value, used in previous studies, 
is somewhat under debate as a harder index can be possible, 
for instance, if the u nderlying population of relativistic elec- 
trons is very narrow ( Katarzviiski et al.l l2006t iTavecchio et alJ 
1200 9*: Lefa et al. 'lOllVh in the case of internal photon absorp- 
tion (Aharonian et al. 2008b ), or in proton-synchr otron models 
(e.g. lAharoni an 2000; Zac haropoulou et ani201 ih . For the lat- 
ter criterion, the expected cascade emission is investigated and, 
additionally, the total intrinsic luminosity is compared to the 
Eddington luminosity of the AGN. Limits on the EBL density 
are derived using each of the criteria individually and for com- 
binations of the criteria. In addition, the influence of individual 
data sets is tested. The obtained constraints reach from 0.4 (im to 
100 ^im and are below 5 nWm"^ sr"' between 8 ^m and 31 [im 
even though more conservative criteria are appUed and the evo- 



lution of the EBL with redshift is accounted for. In the optical, 
the EBL density is limited below 24 nW m"^ sr"' . 

The limits forecast a low level of the EBL density from 
n ear to far in frared waveleng ths also predicted by th e models 
of lKneiske & D ole (201d|) and lDommguez et"al] (1201 Ih which is 
in accordance with MR07. Furthermo re, the co n straint s exclude 
the direct measurements of Matsum oto et a D (l200l . Certain 
mechanisms, however, are discussed in the literature that ef- 
fectively reduce the the attenuation of y-rays due to pair pro- 
duction. For instance, if cosmic rays produced in AGN are not 
deflected strongly in the intergalactic magnetic field they could 
interact with the E BL and form VHE y-ray s that contribute to 
the VHE spectrum dEssev & Kusenkol 120101: lEssev et alj|2010l 
201 li) . Other suggestions are more exotic as the y invoke the con- 
version of photons into axion like particles (e.g lde Angelis et aP 
l2009l: lM irizzi & Montanino 2009) or the violation of Lorentz in- 
variance (e.g. Jacob & Piran 2008). 

Future simultaneous observations of extragalactic blazars 
with the Fenni-LAil and lACTs have the potential to further con- 
strain the EBL density. 
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Fig. 8. Upper limits of this work together with previous limits and EBL models. 
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Appendix A: Evolution of the EBL with redshift 



equat ion for the electron d istribution N(y) in the steady state 
limit (lTavecchioetal.ll201lh . 



N(y) = ^ J dy'Qiy' 

7 



(B.l) 



In this Appendix the influence of the evolution with redshift of 
the EBL density on the upper limits is investigated. For this pur- 
pose, the limits are calculated for four VHE spectra, namely ,^ 
lES 1101-232, H 1426-1-428, MknSOl, and 3C279 which is the Q{y) = (gM^.J) _ 

source with the largest redshift of z = 0.536 in the sample. 
The evolution is in cluded with a phenomenological ansatz (e.g. 
iRaue & MazinI I20b8ft . The four spectra are tested against the 
VHE-HEIndex , PileUp and VHEConcavity criteria and the en- 
velope shape is determined with and without the evolution with 
redshift. Not surprisingly, the EBL density is less confined if the 
evolution is accounted for, as seen from Figure IATTI The differ- 
ences are most pronounced in the optical, where the influence of 
3C 279 is the strongest (light blue shaded region in Figure IaTI i. 
Without taking the evolution into account, the limits are over- 
estimated by up to 40 % at 0.6 \\m. At higher wavelengths, the 
difference is not as distinct as in the optical. This outcome em- 
phasizes that the evolution of the EBL density with redshift has 
a non-negligible effect on upper limits derived from VHE AGN 
spectra, especially for sources with a large redshift. 



where lyl is the energy loss due to IC scattering with CMB pho- 
tons in the Thomson regime with an injection rate 



(B.2) 



The approximation E - Im^c^y is used, with m,, the mass of 
the electron. The cascade spectrum for scattered photons with 
energy e off the e^e~ pair s is then readily calculated to be 
(iBlumenthal & Gouldl[T970l) 



F(6) 



9 e m^c 
64 mcmb(z = 0) 



max[ Ve/4ecMB,rdii.ycng] Zm^c^y 



ncMB(e',z = 0) , 



F'^(e,e\y), 



(B.3) 



Appendix B: Cascade emission formulae 

The electron spectrum due to the interaction of VHE photons 
with the EBL is obtained by solving the corresponding kinetic 



where mcmb(z = 0) = 0.26 eV cm"-', ecMRiz = 0) = 634 (xeV and 
«CMB(f',z) denote the energy density, mean energy and differ- 
ential photon number density of the CMB, respectively. The in- 
verse Compton kernel for scattering on an isotropic photon field 
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Fig. 5. Upper panel: Example of an EBL shape excluded by 
Mkn 501 with the VHEConcavity criterion. Lower panel: The 
spectra of Mkn 501 and Mkn 421 corrected with this particular 
EBL shape. The flux of the latter is scaled by 10"^ for better 
visibility. For Mkn 501, a double broken power law provides the 
best description with a spectral index - -35 at highest ener- 
gies, the maximum value tested in the fitting procedure. In the 
case of Mkn 421, a simple power law suffices. 

in the Thomson regime (Klein-Nishina effects can be neglected 
for the energies of the primary photons considered here) is 



F^^ie, e', y) = 4e'y (2e In e + e + 1 - 2e) , 



(B.4) 



with Q <F^^ <l and e 



e/ (4e'7^)- The lower limit for the inte- 
gration over y is the maximum of three different constraints on y. 
The first one stems from kinematic constraints of Compton scat- 
tering. The second denotes the y-factor for which the electrons 
are deflected outside the opening cone of the blazar jet with an 
opening angle 9c ~ I/Fl, with Fl the bulk Lorentz factor of the 
plasma of the jet. And finally, the third factor gives the minimum 
poss ible Lorentz factor i f the source is active for a certain time 
(see iDermer et al.l 1201 1[ for further details ). However, ydfl and 
yeng are calculated bv IPermer et al.l (1201 1 ^ under the approxi- 
mation of small deflection and observing angles. In the case of 
isotropic emission, i.e. Bigmf ~ lO"'-' G and a lifetime of the 
source of Af > 10^ years, this approximation does not longer 
hold and the lower integration bound is replaced by yfejAecMS- 
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